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, I ^ 

The theory of Singular Lines, constant damping ratio-constant 
undamped natural frequency lines, is derived. 

A limitation of the parameter plane method for characteristic 
polynomials whose coefficients are linear functions of two variable 
parameters which results in undetermined roots is described. The 
addition of singular lines to the parameter plane diagram specifies 
these roots allowing solution for all roots of a given polynomial. 

A general method of solving for singular lines is developed and 
rules for predicting the existence of such lines are stated. A computer 
program which solves for and graphically displays singular lines in 
addition to constant zeta, omega, and sigma loci is presented. 

Singular lines are considered in terms of dominance and macro- 
scopic root sensitivity. A dominant root line in the parameter plane 
is illustrated. 

Examples which demonstrate the application of singular line 
theory to linear, multivariable control systems are presented. 
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INTRODUCTION 



Prior to 1959 analysis and synthesis of control systems were 
carried out primarily by root locus and frequency response techniques. [1] 
Mitrovic*s Method [2], which was introduced in 1959, allowed the deter- 
mination of a system’s root locations as two coefficients- of the 
characteristic equation were varied. The method specified the variable 

coefficients, and B as the two lowest order terms of the character- 
1 0 

istic equation. For example: 

a s^ + a ^ + + a^s^ + B-s + B =0 (1-1) 

n n-1 2 1 0 

where: 

a, k = 0, 1, . . . ,n are real constants 

k 

B^ and B are real variables 
1 0 

Application of the method resulted in constant zeta, omega, and sigma 

curves plotted on the B^ - B^ plane specifying the roots of equation 

1-1 for any choice of C and o) . The Coefficient Plane Method includes 

n 

Mitrovic’s original work and extension of his method to include variation 
of any two of the coefficients of the system characteristic equation. 

In 1964, Siljak [3] extended Mitrovic’s work to control systems 
in which two adjustable parameters appeared linearly in the coefficients 
of the characteristic equation. In the linear system case, Siljak 
considered the characteristic equation: 

n , 

f(s) = I as = 0 (1-2) 

k=0 
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where : 



“k ■ \ 



(1-3) 



and: b, , c, , d, are real constants 

k k’ k 

a, 3 are real, variable parameters 

This method produced constant zeta, omega, and sigma curves plotted on 
the a 3 plane specifying the roots of equation (1-2) for any choice of 
a and 3. The Parameter Plane Method, as Slljak named it, is a simple 
procedure for factoring characteristic polynomials and displaying the 
results in a parameter plane diagram. 

The Parameter Plane Method was extended In 1965 to a general case 
in which coefficients of the characteristic equation are a nonlinear 
combination of two adjustable system parameters. [4], [5], [6] Specif- 
ically, the product case where the coefficients, a^^, of equation (1-2) 
were: 

aj^ » bj^a + Cj^3 + h^a3 + d^ (1-4) 

where : 

b, ,c, , h, and d, are real constants 

k k’ k k 

a, 3 are real variable parameters 

A further extension of the method to include characteristic equation 
coefficients of quadratic form is limited to the specific case of third 
order systems. [7] 

The Parameter Plane Method was conceived as an approach to the 
analysis and synthesis of feedback control systems. It allows the 
designer to obtain information about system relative stability and the 
effect of parameter adjustments on stability. By adjusting pole-zero 
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locations in the system transfer function the designer maintains control 
over both transient and frequency responses. [8] The method is 
particularly applicable to multiparameter, multiloop control systems 
with more than one adjustable parameter since it is based on the study 
of the system characteristic equation written in a general form. In 
general the parameter plane method may be applied effectively to any 
engineering problem in which it is necessary to determine how variations 
of parameters in the characteristic equation effect root locations. 
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II 



THE PARAMETER PLANE METHOD 



This section will review Siljak’s derivation of the parameter 
plane method for the linear case [3], discuss the use of the parameter 
plane for system stability analysis, and illustrate the primary rules 
for mapping points and curves from the s-plane to the a g plane. 

Parameter Plane Equations 

Consider the characteristic equation: 

/ 

f(s) = I = 0 (2-1) 

k=o 

where: 

the coefficients a^ (k=o,l, n) are real and s is 

the complex frequency variable: 

s = a + jo) = "C (jo + jo3 (2-2) 

n n I 

where to is the undamped natural frequency and c; is the damping 
n 

ratio 

By forming powers of s it can be shown [3] that: 
where: 

\(-C) = (-1)*" (2-4) 

Uj ^(-0 + (- 1)^‘^\(0 

The Uj^ and Tj^ are Chebyshev functions of the first and second 
kind respectively. They are given by the recursion formulae [3]: 
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- 0 
- 0 

where: 



(2-5) 



Tq(C) = 1 Uq(;) = 0 

= C = -1 



Substituting equation (2-3) in equation (2-2) and equating the real and 



imaginary parts of the resulting equation to zero independently gives : 



I 



k=0 



“k"„ 



= 0 



( 2 - 6 ) 



t “k“n “ ''k<-«> - “ 



From equations (2-5) we obtain the interrelation equation: 

Substitution of equations (2-7) and (2-4) in equations (2-6) produces, 
after simplification, equations in one Chebyshev function: 



I “k“„ " “k-i<^) = “ 



(2-8) 



I v„ " “ 



Assume, as Siljak did, that the coefficients a^ of the character- 
istic equations are linear combinations of parameters a and 6 of the form: 



\ ■ *’k“ 



+ c, 6 
k 



+ d. 



(2-9) 
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where : 



b. , c, and d, are real constants 
k k k 



a and 3 are variable 



Substitution of equation (2-9) in equation (2-8) yields: 



Bj^a + Cj^3 = 

B^a + 0^3 = -^2 



(2-10) 



where : 



k=0 



<=2 ■ X < 2 - 11 ) 

k=0 k=0 



n 



n 



“l - I <-«"\“n\-l<') “2 - I 

k=0 k=0 



The functional dependence of B^^, B^i C^, D^, and on ^ and o)^ is 

omitted in equations (2-10) and succeeding equations for simplicity of 
notation, A Table of the Chebyshev functions of the second kind, 
is given in Appendix I, The numerical value of U^(^) selected C is 
given in Appendix II. 

Application of Cramer’s Rule for solution of simultaneous linear 
equations to equation (2-10) yields the desired parameter plane solution 



equations : 



a 



^l°2 ~ 

®1^2 ■ ®2^1 



6 = 



^ 2^1 ~ h ^2 
- «2^1 



(2-12) 
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As defined by Siljak [3] the parameter plane or a3 plane is a 
rectangular coordinate plot with a the abscissa and 3 the ordinate. 



Equations (2-12) give a and 3 as functions of C and Fixing ^ ^ 

and varying from 0 to infinity in equations (2-12) produces a curve 
in the a3 plane, a constant zeta curve. This curve specifies the a3 
pairs which will cause equation (2-1) to have a pair of complex roots 
with the required damping ratio Similarly, fixing and varying 

C from -1 to +1 in equations (2-12) produces a curve in the parameter 
plane, a constant omega curve, specifying the a3 pairs which will cause 
equations (2-1) to have a pair of complex roots with the required 
natural frequency Equations (2-12) thus allow mappings of points, 

excepting real axis points, from the s-plane to the parameter plane. 

For real axis points in the s-plane, replace s in equation (2-1) 
by s = -a to obtain: 



n , 

f(s) = 1 a (-a) = 0 

k=0 



(2-13) 



Now, substitute equation (2-9) for to obtain: 



n 



\ (bj^a + + dj^) (-0) =0 



k=0 



(2-14) 



Simplifying equation (2-14) yields: 
aB (a) + 3C (a) + D (a) = 0 
where: 

B(0) = I (-I)'" b.a^ 
k=0 



(2-15) 



r k k 

C(a) = I (-1)*" d o*" 

k=0 



(2-16) 



n 

I 

k=0 



D(o) = I (-1)^ d 
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For a given value of a the functions B(a), C(a), and D(a) are constants 
and equation (2-15) is the equation of a straight line on the a3 plane, 
the locus of a3 points corresponding to real roots s = -a, 

I • 

Stability Analysis and Mapping 

Absolute stability analysis of a linear control system consists 
of determining the existence of roots of the characteristic equation in 
the right half of the s-plane. Similarly, relative stability analysis 
consists of determining the existence of roots within specified areas of 
the s-plane, generally defined 'by constant C or constant o)^ loci. 

Mapping of the specified areas from the s-plane to the a3 plane using 
equations (2-12) and (2-15) allows a designer to choose or adjust para- 
meter values so that characteristic equation roots lie within required 
areas thereby ensuring stability requirements are met. The rules and 
graphical techniques for the mapping process and the properties of the 
mapping are discussed in detail by Siljak [3] and Thaler [9]. The inter- 
pretations of the curve on the parameter plane and determination of 
stability is not a simple process; the reader is again referred to 
references [3] and [9] for detailed information. 

The major mapping rules will be illustrated in the following 
example. 

Example I ; 

Consider a linear system whose characteristic equation is the 
fourth order polynomial: 



f(s) = 0.04 + 0.34 + (0.2a + 1.12)s^ + 

(0.5a + e + 1.7)s + 2e + 1 = 0 



(2-17) 
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From the parameter plane diagram for the polynomial shown in Figure 1 
we can make the following observations: 

1) An s-plane complex point maps into a single a3 plane point. 

2) An s-plane real axis point maps into a straight line, a 
constant sigma line, on the ag plane. Real roots may be 
evaluated from the a lines. For example the point M(2. 0,0.83) 
on the a = 1.5 line has a root s = -1.5. 

3) A constant damping ratio, natural frequency or settling time 
contour maps into a constant r,o) or ^o) curve on the ag 
plane. Complex roots of the polynomial for specific values 
of a and g are determined by the ^ and curves. For 
example point M(2.0, 0.83) is at the intersection of the 

C = 0.5 and = 3.0 curves, therefore the polynomial has 
complex roots at s = -(0.5) (3.0) + j 3.0~^l-(0.5)^ 

= 1.5 + j 2.6 

4) Constant sigma lines on the ag plane are tangent to the 

C = 1 curve at 0 ) =1.5. 

n 

5) An s-plane stable area maps into an area of the ag plane in 
which all the roots of the polynomial have negative real 
parts, i.e., into a stable root area. For example, consider 
the area in the left half s-plane bounded by the origin, the 
radial ? = 0.5 lines, and the infinity of the s-plane. 

Setting s = 0 in equation (2.17) yields g = -0.5 as the a = 0 
boundary in the ag plane. The C = 0.5 radial lines map into 
the ? = 0.5 contour of Figure 1. The enclosed area on the 
plane of Figure 1 is, therefore, the desired stable root area. 
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PARAMETER PLANE EXAMPLE 
^TH. ORDER POLYNOMIAL 
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BASIS OF THIS STUDY 



This section introduces a heretofore unknown limitation of the 
parameter plane method which results in undetermined roots, illustrates 
the limitation with a specific example, describes initial investigation, 
and then formulates a hypothesis as to the cause of undetermined roots. 

Sections IV and V develop the mathematical theory, solution 
method, and computer implementation extending the parameter plane method 
to ensure determination of all roots of a given polynomial. Sections VI 
and VII consider this extension of theory in terms of dominance and root 
sensitivity. Section VIII contains examples of analysis of specific 
linear control systems which indicate the applicability and potential 
of the theory developed. Section IX concludes this study with comments 
on work completed and recommendations for further investigation of para- 
meter plane theory and methods. 

Limitation of Parameter Plane Theory — Undetermined Roots 

The parameter plane method derived by Siljak [3] and the computer 
program, PARAM A, written by R. M. Nutting [8] purportedly solved the 
problem of determining all roots of a given polynomial, characteristic 
equation, in terms of variable parameters a and 6 displaying the results 
on the parameter plane as constant ?, o)^, and curves and constant a 
lines. In most cases the existing parameter plane method does, in fact, 
solve for all roots of a given polynomial, and in all cases such infor- 
mation as it does produce is correct. There are, however, situations 
in which the present parameter plane method does not provide a solution 
for all roots of a polynomial nor correctly predict the existence of 
complex roots in specific areas of the a3 plane. Example II is an 
illustration of such a situation. 
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Example II 



Consider a linear control system whose characteristic equation is 

+ 80s^ + (20a + 1600)s^ + 840as^ + 

2 (3-1) 

(1600a + 4003) s + 16003s + 16003 = 0 

Solution by the parameter plane method results in the parameter plane 
diagram shown in Figures 2 and 3, The use of two figures was necessi- 
tated by the complexity of the constant C curves. Although it is not 
obvious in Figures 2 and 3 the constant ^ curves all go to the infinity 
of the parameter plane as increases. Although the ^ = 1 curve has 
discontinuities it effectively circles the origin clockwise and then goes 
to infinity in the second quadrant of the ag plane. The computed values 
of a and 3 for the constant ^ curves and the use of smaller graph scale 
for 3 confirms these statements. 

The area of absolute stability, all roots in the left half s-plane, 
is bounded in the a3 plane by the C = 0 curve, the positive a axis and 
infinity. 

In this example the present parameter plane theory correctly 
solves for all roots of the polynomial in the complex root area of the 
a3 plane in Figure 2. In the real root area bounded by the C = 1 curve 
and the stable area limits the present parameter plane method does not * 
determine all roots of the polynomial; only two or four real roots are 
determined by constant o line intersections and no information is 
presented for the remaining roots. For example, for point (40, 300) 
of Figure 2 the roots are 
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= - 2.12 



s., = -37.88 s, = ? 

2 4 



Similarly, for point M 2 (20, 26.6) the roots are 
= - 2.77 
S 2 = - 9.77 
83 = -27.5 
= -39.3 




Existence of Undetermined Roots in Parameter Plane Solutions 

Further investigation of the situation outlined in Example II 
yielded the following results. In that part of the parameter plane for 
which present theory predicted only real root^ specific a - 6 points 
were found to have either one or two pairs of complex roots. In 
addition, certain points had the same complex pair in common. For 
example, for point (40, 300) of Figure 2 the complex roots are: 

~ “ 1*06 - j 1.83 
= -18.94 - j 31.24 

For points (20, 26.6) and (80, 113) the complex roots are 

s^^g = - 0.35 - j 1.16 

The results were obtained by substituting specific a3 pairs in equation 
(3-1) and solving for the roots of the resulting equation. Moreover, 
joining points M 2 and M^ by a straight line and investigating specific 
points along the line showed that all such points had a common pair of 
complex roots. Thus a line of constant ^-constant o)^ was located on the 
parameter plane. 
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As previously stated, the existing parameter plane method correctly 
solves for all roots of a given polynomial in the majority of cases and in 
all cases the root values which it does determine are correct. It is 
therefore obvious that existing theory as derived [3] is correct insofar 
as it goes, but it overlooks or dismisses some special case. Reviewing 
the derivation, Cf. ante pp. 12-16, it was apparent that the use of 
Cramer's rule in solving equations (2-10) presumed that the coefficient 
matrix was non-singular, i.e., that the value of the determinant which 
formed the denominator in relations (2-12) was non zero, for all cases 
of interest. This then seemed to be the part of existing theory which 
could lead to undetermined roots, i.e., incomplete solution of the 
problem, and this hypothesis formed the basis for the remainder of this 
study. 
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IV 



SINGULAR SOLUTIONS OF LINEAR 
PARAMETER PLANE PROBLEMS 



In Section III a heretofore unknown limitation of the parameter 

plane method resulting in undetermined roots was described. Example II 

presented a specific illustration of the problem of undetermined roots. 

It was shown that at least one constant C - constant (o line existed and 

n 

that it provided solutions for some of the roots which existing para- 
meter plane theory could not determine. It was then hypothesized that 
the assumption that use of Cramer’s Rule for solving the simultaneous 
linear parameter plane equations, equations (2-10), produced the complete 
solution for all cases of interest excepting that of real roots was in- 
correct. This section extends parameter plane theory to include the 
complete solution of equations (2-10) by solving for presently undeter- 
mined roots. The theory of constant C *“ constant and lines defined 
SINGULAR LINES is derived, a tractable singular line solution method is 
developed, and general rules for determining the existence of such lines 
are stated. 



I. SINGULAR LINE THEORY 



Mathematical Basis 

The solution of a system of linear non-homogeneous equations by 
matrix methods is well known and is contained in any standard textbook 
of linear algebra [10] or matrix textbook [11]. Consider the system of 
equations : 





h2 •“ 


3ln 




1 

1 — 1 
X 

1 




i 

1 


^21 


322 * • • 


32n 




^2 


= 


h2 




3m2 


a 

mn 




X 

n 




h 

n 
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or, in compact notation, 



AX = H 



(4-2) 



By way of review the following definitions and statements are given 

without amplification or proof: 

1. A = [a,.] is the coefficient matrix. 

ij 

2. [AH] is the augmented matrix. 

3. A system which has a solution is said to be consistent. A consistent 
system has just one solution, a unique solution, or infinitely many 
solutions . 

4. A system which has no solution is said to be inconsistent. The 
problem leading to such a system of equations is indeterminate. 

5. In the consistent case, the coefficient matrix A and the augmented 
matrix [AH] have the same rank. In the inconsistent case, they have 
different ranks. 

6. When one or more equations of a system can be derived from another 
by multiplication of all of their terms by a constant, the equations 
are dependent or equivalent. 

7. The matrix A is called non-singular if its rank r = n, that is, if 
|a| ^ 0. Otherwise, A is called singular. 

8. The system of equations has a unique solution provided the common 
rank of matrices A and [AH] is equal to n, the number of variables, 
that is provided |a| ^ 0. Thus the system has a unique solution if 
A and [AH] have the same rank and if A is non-singular. 

9. In a consistent system of rank r < n, a solution can be obtained for 
r variables in terras of the remaining n - r variables. 



26 



Derivation of the Singular Line Solution Method 

Consider the linear parameter plane equations, equations (2-10), 
which comprise a system of linear non-homogeneous equations. Transposing 
terms and rewriting the equations gives 



B^a + C^3 = - 



(4-3) 



The coefficient and augmented matrices for equations (4-3) are 



A = 



AH = 




(4-4) 



(4-5) 



The present parameter plane method solves equations (4-3) for those 
cases in which the common rank of A and [AH] is two as follows: 



a - 



3 = 



I— 1 

Q 

1 




.‘°2 




1 — 1 
PQ 




_«2 








_«2 


■'' 2 _ 


1 — ! 
PQ 




_®2 





‘^ 1^2 ^ 2^1 



- ® 2^1 



® 2°1 - " l °2 



- ®2"l 



(4-6) 
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There remain, however, the cases in which equations (4-3) are 
linearly dependent. In these cases matrix A is singular rank one; 
therefore, the above method of solution is not applicable since the 
determinant of A, |a| , equals zero. To obtain a solution we must first 
choose a C " pair such that matrices A and [AH] have common rank one 
and then solve one of equations (4-3) for one unknown, 3, in terms of 
the other, a. 

For matrices A and [AH] to have common rank one of the following 
conditions must apply; 




V2 - “ 



<=l“2 - ■ 0 



“2“l - ®l‘>2 ■ “ 



( 4 - 7 ) 



( 4 - 8 ) 



( 4 - 9 ) 



To choose a C 
determinant 



- 0 )^ pair which satisfies these conditions, expand the 



®2 S 



in terms of ? and o)^, substitute a specific 



value of C = equate to zero, and solve the resulting polynomial for 

0 )^. Real roots of the polynomial and the specific value of C = are 

then substituted in equations (4-8) and (4-9) ; those values of which 

satisfy the equations, o) , are the values for which matrices A and [AH] 

s 

are both rank one. 

Substitution of a C - o) pair in equations (4-3) and solution 

s 

of the resulting linearly dependent equations for parameter 3 in terms 
of parameter a yields 
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(4-10) 




where: 



k=0 s 

k=0 s 

'>1 ' 

k=0 s 



“2 ' I 

k=0 s 

k=0 s 

“2 - J„ <-» Vn 

k=0 s 

(4-11) 



Equation (4-10) is the desired singular line solution. It is the equa- 
tion of a straight line on the ag plane with slope - and g axis 

intercept - points on the line are maps of the s-plane point 

(C i m ) thus it is a constant L - constant ca line. Substitution of 
s n n 

s 

any values of a and g corresponding to the infinity of a-g pairs defined 
by equation (4-10) in the characteristic equation (2-1) results in 
s-plane roots: 




The derivation leading to equation (4-10) is correct, but it is 

by no means a practical or appealing solution method since the procedure 

for solving for ^-cd pairs which result in singular lines (c , m ) is 

n s n 

s 

both tedious and time consuming. The presence of Chebyshev functions 
in equations (2-11) produces polynomials of order n and n-1 in where 
n is the order of the given characteristic equation (2-1). Expansion 
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of equations (4-7), (4-8) and* (4-9), which in itself is a tedious 
process, leads to a polynomial of order 2n-l in For complex systems 

with higher order characteristic equations the resulting polynomial may 
be extremely difficult to solve. 

In order to derive a more tractable solution method consider a 
characteristic equation of the form: 

n , 

f(s) = I ih.a + c, 6 + d, )s‘^ = 0 (4-12) 

k=0 ^ 

Setting n=i for i=2,3; computing B^, C^, and in equations (2-11); 

expanding equation (4-7); and grouping terms yields: 

For n = 2 : 



For n = 3: 



V2-®2^1 = 



[bQCi-biCojoJn - [bQC2~b2CQ] 2?(o^ + 



= 0 



Details of the expansion for n=2, 3, 4 and the extension to the general 
case are contained in Appendix III. The result for the general case is 
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B. 



Wr 






(4-13) 



or, in compact notation: 



® l‘=2 - 



Vl ' I. 



ce . . 
iJ 



(4-14) 



where: 



CE = [ce^^] is an (nxn) upper triangular 
matrix whose elements are 

[b. iC ^c. ] (-1)^ ^U.(c)w 3 2 ^ ^^2 

1-1 i+j-1 i+j-1 1-1 j n 

(4-15) 

0 i+j ^ n+2 

Equation (4-14) results in a polynomial in C and o)^ of order 
2(n-l) in Substitution of a specific C = and solution for real 

values of gives the s-plane coordinates for which the coefficient 
matrix is singular, equation (4-4). The real values of o)^ and the 
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specific ? * which satisfy equations (4-8) and (4-9) are the 

coordinates (t , 03 ) that ensure matrix A is singular and matrices A 

s 

and [AH] have rank one. ? and o) are substituted in equation (4-10) 

s 

which is now a linear equation in two unknowns a and 3. The resulting 
equation defines a line of constant C " constant a singular line, 
on the a3 plane. 

The above solution method appears cumbersome when described but, 
in fact, it proved to be very simple to apply. To illustrate the solu- || 

tion method Examples III and IV are given below. In common with the 
present parameter plane method, the complexity of computation increases 
with the order of the given polynomial and computer solution is preferable 
for higher order-systems. 

Example III 

Consider the characteristic equation given in Example I: 

f(s) - 0.04s^ + 0.34s^ + (0.2a + 1.12)s^ + (0.5a + g + 1.7)s + 28 + 1 - 0 '• 

(4-16) • 

• i 

H 

I 

For ease of coefficient component identification, display the terms as i' 

follows : ■ 




Apply equation (4-14) to compute: 
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Simplifying and equating to zero: 



0) ^ - 4^0) +5 = 0 

n n 



Choosing C =0.5, for example, and solving for the values of o) for 
s n 

which matrix A is singular: 

2 

0) - 2^0) +5 = 0 

n n 



= 2 - A- 20 
2 




Since there are no real values of the matrix is non-singular for all 
values of no constant C-constant o)^ line exists, and the present 
parameter plane method correctly solved for all roots of the given poly- 
nomial. 

Example IV 

Consider the characteristic equation given in Example II: 
f(s) = + 80s^ + (20a + 1600)s^ + 840as^ + (1600a + 400g)s^ + 

16006+ 1600;= 0 (4-17) 

Display the terms as: 
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Apply equation (4-14) to compute: 






- Co] 

-Cl600^]“Af 

-[+oo(840jq^ 



+ [lt»o(840)]2$o;+ -[|«OoM](4S*-0«^* 

+ [400 (xo)] zt 



Simplifying and equating to zero: 

^ - (8c^+19)oj ^ + (16e^+160?)u) ^ - (336c^+76)u) + 320? = 0 

n n n n 

Choosing ? =0.5, for example, and solving by synthetic division for u) 

s n 

gives: 

0) ^ - 42to ^ + 1640) ^ - 320o) + 320 = 0 
n n n n 

0 )^ = 2.1115, 37.889, and a complex pair 

Thus the matrix A is singular for ? = 0.5, = 2.1115 and for ? =;0.5, 

u) = 37.889. 
n 

Substituting the first pair in equations (4-8) and (4-9) leads to the 

results 0=0 confirming that (? ,m ) = 0.5, 2.1115 is a singular point. 
From equations (4-11) : ® 

= (-l)^(1600)(2.1115)^(l)+(-l)^(840)(2.1115)^[2(0.5)]+ 

(-1)^(20) (2.1115)^[4(0.25)-1] 

= -774.06 



34 



= (-1)^(1600) (2.1115)*^(-l)+(-l) (1600) (2.1115) (0)+ 

(-1)^(400) (2.1115)^(1) 

= M83.4 

= (-1)^(1600) (2. 1115)^ [4(0. 25)-l]+(-l) ^(80) (2. 1115) ^[8(0. 125)-2]+ 

(-1) ^(2. 1115) ^[16(.5)^-12(. 5)^+1] 

= 3268.8 

Substituting these values in equation (4-10) yields 

^ 3268.8 774.06 

^ 183.4 183.4 “ 

= 17.81 + 4.223 a 

which is a line of constant ^“constant o) (.5, 2, 1115) whose g axis 

n 

intercept is 17.81 and whose slope is 4.22. Similarly, the values 
^ = 0.5, 0 )^ = 37.889 satisfy equations (4-8) and (4-9) and define a 
singular line: 

6 = 5740 + 75.777a 

Mitrovic^s Method and the Coefficient Plane Method 

Singular line theory must be compatible not only with the para- 
meter plane method but also with its special cases, Mitrovic’s method 
and its extension to the coefficient plane. 

The coefficient plane method applies when any two coefficients of 
the characteristic polynomial, equation (2-1), are variables and 
where n ^ p > q ^ 0. The solution parallels that for the parameter 
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plane, outlined in section II, up to and including equations (2-8). 

Substitution of variable coefficients A and A in equations (2-8) yields 

p q 

(-1)»A yu (?)+(-1 )’a U.’U (C) - - I 

p.q 



n 



(-l)PAp(oPUp(O + (-l)'lA^a)%^(c) = I (- 1 )\j^o)NJj ^(0 

lc~0 

^ p.q 



(4-18) 



Considering these as a system of linear, non-homogeneous equa- 
tions in two unknowns, A and A , the coefficient matrix is 

P q 



(-l)Po)Pu^ . (?) (-I)VU ,(?) 

p-1 q-1 



(-l)Po)PUp(?) 



which simplified is: 

(-w)"‘i 



(-i)Vu^(?) 



(4-19) 



(-(j)P”% 



P-1 ^-1 



(-uj)P-'^U 



U 



(4-20) 



Considering that p > q, inspection of a representative table of Chebyshev 
functions of the second kind. Appendix II, shows that the determinant of 
the coefficient matrix, equation (4-9), cannot equal zero. Thus, since 
the matrix is non-singular except for the trivial case cd^ = 0, equations 
(4-18) cannot be linearly dependent and the solution for A^ and A^, if 
it exists, must be unique. 

In order to test the applicability of singular line theory to 
coefficient plane problems, recast them as parameter plane problems. 

All possible cases are covered if we consider: 
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A (or A ) = f(a) 

p q 

(4-21) 

A (or A ) = f(6) or f(a,3) 

q P 



In any case, since only two coefficients of the characteristic equation, 

(2-1), contain the variable parameters the expansion of the coefficient 

matrix, equation (4-13), contains one and only one term. The result is 

that singular line theory always produces the solution for the trivial 

case, (A) =0, and no singular lines are defined. Singular line theory 

s 

is, therefore, compatible with the coefficient plane and Mitrovic’s 
Method, 



II, CONDITIONS FOR EXISTENCE OF SINGULAR LINES 

In order for singular lines to exist the linear, non-homogeneous, 
parameter plane solution equations, (4-3), must be dependent. The 
necessary conditions for this are that their coefficient matrix must be 
singular and their coefficient and augmented matrices must have common 
rank one. 

Ideally an existence theorem which clearly defines the conditions 
on the characteristic polynomial necessary for singular lines should be 
stated. This cannot, however, be done since the coefficients and 
constants in the parameter plane solution equations, equation (4-3), are 
each involved functions of n+2 quantities where n is the order of the 
characteristic polynomial, that is, functions of b^^ or c^ or d^ and C 
and u)^. Certain general rules which indicate the likelihood of the 
existence of singular lines can, however, be formulated. They are 
1. The variable parameters a and 3 must be contained in more than two 
of the coefficients of the system characteristic equation. 
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2. Equation (4-13) indicates that for a characteristic equation of 
order n: 

a. For each value of C there is a maximum of 2(n-l) values of co 

n 

which produce singularities, excluding the trivial case of o)^=0, 

b. For each value of there is a maximum of (n-1) values of C 
which produce singularities. 

There are an infinity of possible values of C for 0 < C £ 1, each 
associated with a maximum of 2 (n-1) values of o)^ which produce 
singularities. Thus, if the conditions for singular line existence 
are met, there is a continuum of such lines spanning the real root 
section of the parameter plane stable area. 

3. Coefficients of a and 3 in the higher powers of s in the characteristic 
equation have the most effect on the possible existence of singular 
points. 

4. If in the characteristic polynomial ^ the 

maximum possible number of values of per constant C curve which 
produce singular points is M-1 excluding the trivial case of o)^ = 0. 

5. If in the parameter plane diagram the stable area is clearly split 
into a complex root area in which all the constant C contours 
eventually go to infinity and a real root area which contains only 
constant a lines, the possibility of existence of singular lines 
is high. 
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V 



COMPUTER PROGRAMMING OF THE SINGULAR LINE SOLUTION METHOD 



Appendix IV is a listing of the digital computer program developed 
in the course of this study. The program, PARAM S, is comprised of the 
major parts of PARAM A, a program for solution of the linear parameter 
plane problem written by R. M. NUTTING [4] , and a new sub-program to 
implement the singular line solution method developed in section IV of 
this study. Computation and graphing of constant C curves, constant o)^ 
curves, constant a lines and constant ^ - constant o)^ lines, singular 
lines, is provided. The many options of quantities computed, curves 
plotted, and output data printed and the format of necessary input data 
are explained in the comment sections of Appendix IV. ' 

The computation procedure used in the singular line sub-program 
for a characteristic polynomial of the form: 



f(s) = I a.s = 0 

k=0 



(5-1) 



where: 






and for a specified value of damping ratio, Q, is as follows: 

1. Read input data and the coefficients of the characteristic poly- 
nomial where k = 1, 2, n+1. 

2. Form the matrix CE = [ce..l . , ^ 

ij" 1=1, 2,.,., n 

j — 1» 2 , « e « ,n 

whose elements are 



(b.c -b c.) 

1 i+j i+j 1 



i+j < n+2 
i+j ^ n+2 



39 



3. Form the matrix CF = 

. . . ,n 

whose elements are 

1 =■ j 

0 i ^ i 

CF is a signed diagonal matrix whose elements are formed from the 
Chebyshev functions of the second kind for the chosen value of 5 . 

4. Generate the AA matrix where: 

AA = CE X CF 

Sum the elements of the AA matrix. The resulting sum is the expan- 
sion of (B^C 2 ~B 2 C^) , equation (4-13). The summation process groups 
the terms in powers of m and stores the coefficients of w In 
ascending order. 

5. Solve for using the quadratic formula or subroutine POLYRT as 

appropriate. POLYRT is a FORTRAN 60 subroutine capable of solving 

for the roots of a polynomial, order n where 2 < n < 100 , with 

arbitrary complex coefficients. The subroutine uses Lehmer's 

method [ 12 ] to approximate a root and then improves on it with 

Newton's method. POLYRT is a relatively slow but accurate subroutine 

which requires no supervision or Initial guesses of root values. For 

a polynomial with repeated real roots, the results often have 

-3 

erroneous complex parts which are generally less than 2 . 0 x 10 

6 . Test the output of POLYRT retaining only real roots; these, together 
with the specified value of c, define the singular points. If there 
are no real roots print "Matrix Non Singular", and halt. 
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?• For each real value of o)^ and the specified C generate 

and D^, equations (2-11); test that the conditions for 
dependent solution, equations (4-7), (4-8), and (4-9) are satisfied. 

If print "System of Equations indeterminate. Rank 

Zero." re loop for next value of or halt as appropriate. If either 
(4-8) or (4-9) is not satisfied, print "System of Equations Incon- 
sistent. No Solution." re loop for next value of w or halt as 
appropriate. 

8. If conditions are satisfied computer g for values of a determined by 
X scale and lYRIGHT input data, equation (4-10), and plot a graph of 
the singular line. As determined by program options print: 

XAXIS intercept, slope, third parameter, and 
b. roots of the characteristic polynomial at selected points on 
the singular line. 

Subroutine POLYRT is again used in solving for the roots. The printout 
of roots serves both as a check on the correctness of the plotted 
singular line and as an indication of the movement of other roots as 
the singular line is traversed. 

Example V 

Consider the characteristic polynomial: 

+ (20a+42)s^ + (20a+400B+l6l)s^ + (2440a+l600g+3280)s + 

(2360a+l600p+l600) = 0 (5-2) 

The output graph from program PARAM S, the parameter plane diagram, for 

this polynomial is shown in Figure 4. It is comprised of constant ? 

curves, constant curves, and singular lines, constant constant 

lines, for selected values of C and m . 

n 



41 



The computed output data associated with each singular line is 



contained in Appendix V. The two forms of computer printout of singular 
line information provided by program PARAM S are illustrated, Table 111 
contains the values of XAXIS Intercept, and slope for each singular 

line. Table IV Includes the roots of the given polynomial at selected 
points on a representative singular line. 

The constant c curves are computed for values of from ,02 to 
200, Although it is not obvious in Figure 4, computer data printouts 
and smaller scale graphs show that the constant C contours describe 
short arcs in the second quadrant of the plane, jump discontinuously 
to the fourth quadrant, then follow smooth curves through the fourth 
and third quadrants prior to proceeding to infinity as shown, 

The singular lines in the stable area, effectively bounded by 
the ? * 0 contour and the postilve a axis, each define a line of 
constant ^ - constant cd and are coincident with constant o lines, For 



example, singular line 43 for which C *0,5, u) *12.723 is a line of 

B XI 

8 

pairs for which the given polynomial has roots: 



- -6.3614 i jll.OO 
and s ■ -0 

- -0.8077 

As shown in Figure 4 the remaining real root moves from 8*71,5 to s*-28.1 
as a is increased from -5,0 to 0 and the root continues to move to the 
left in the s-plane as a is increased further. The remaining singular 
lines are contained In the unstable region of the plane; that is, one 



n 
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or more roots of the given polynomial has a positive real part for any 
a-6 pair in this region. As indicated on the parameter plane diagram, 
these singular lines are grouped in three distinct sets. Two of the 
three sets contain singular lines in the ag plane for which a specific 
pair of complex roots, the singular pair, are in the right half of the 
s-plane. The other set contains singular lines which specify the 
singular complex pair and one real root in the left half of the s-plane. 
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I 




K'XRLf. - 5.00E-f08 U^aT^/IMCH. 

7‘iCflLE • J.00E-f0i UNIT3^iNCH. 

FIGURE 4. 



SIMGULflR LINES ON THE PARAMETER PLANE 
OF A 4TH. ORDER POLYNOMIAL 
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VI 



DOMINANCE OF SINGULAR ROOTS 



The concept of dominant mode analysis and design has wide applica- 
tion for linear systems. Most system criteria in the s-plane such as 
damping ratio, natural undamped frequency, settling time, and peak over- 
shoot have meaning only for a dominant mode system. The no zero, two 
complex pole system model is accepted as the standard model for s-plane 
analysis of linear systems. The model characteristic equation is 



the characteristic equation approximates equation (6-1) by choosing one 
pair of complex roots to satisfy equation (6-2) and forcing the 
remaining roots to be far to the left in the s-plane or placing zeros to 
cancel undesired poles. The chosen roots then dominate the transient 
response. Since the resulting response is essentially that of a second 
order system, C and o)^ may be taken directly from the dominant root 
location. The peak overshoot, settling time, bandwidth, etc., may then 
be read directly from a set of standard second order system graphs. In 
addition, a third order system model with no zero, two complex poles and 
a real pole in the left half s-plane near the origin may be considered 
to approximate the standard second order model for analysis purposes. 

For systems with higher order characteristic equations the 
approximate analytical method described in reference [13] may be used for 
dominant mode design. Essentially, the approximate method consists of 




(6-1) 



for which the poles are located at: 




( 6 - 2 ) 



When possible a designer attempts to build a system such that 
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writing the system characteristic equation in terms of its parameters 
and factoring the characteristic equation to derive two equations, one 
second order in s. The values of system parameters are then chosen to 
ensure that the roots of the second order equation are dominant, and 
the roots of the other equation are far to the left in the complex 
plane. 

Having designed a dominant mode system, if the values of ? and 

pertaining to a singular line on the parameter plane diagram (c , u> ) 

s n 

+ I 2 ® 

define s-plane roots s ■ -5 uj - Jw which are equal to or near 

"s * 

the chosen dominant roots, then there is a point and possibly many points 
on the singular line which define values of parameters ot and 3 which en- 
sure dominance. For example, singular line 43 of Figure 4, Example V, 

for which 4 ■ 0.5 and co • 12.723, is a dominant root line. Table IV 

s n ’ 

s 

in Appendix V shows that any point on the section of this line in the 
stable area defines values of parameters a and 3 for which the character- 
istic equation (5-2) has the desired dominant roots, a real root near 
the origin, and a real root far removed in the left half of the s-plane. 

When the fortuitous coincidence of chosen dominant roots and 
singular line fixed roots occurs, as in Example V for singular line 43, 
a dominant root area may be shown in the a3 plane. The boundaries of 
this dominant root area, which are determined by the degree of dominance 
required, are divergent singular lines. For example singular lines 33 
and 53 of Figure 4 bound a dominant root area for which domlnancy is 
between 9.5 and 11.0. 
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VII 



SENSITIVITY ANALYSIS OF SINGULAR LINES 



There have been many sensitivity functions defined for the effects 
of parameter variations on the transfer functions of linear control 
systems. The root sensitivity functions developed by Kokotovic and 
Siljak [14] provide the most general solution to the sensitivity problem 
for small parameter variations in linear control systems, and are most 
applicable to parameter plane methods. In addition the macroscopic root 
sensitivity functions defined by F. H. Hollister [5] are applicable for 
large parameter variations. 



Sensitivity Equations 

This formulation is an abbreviated composite version of the 
derivations given by Kokotovic and Siljak [14] and F. H. Hollister [5] , 
and it presents only the salient points. Consider the characteristic 
polynomial: 

I a = 0 (7-1) 

k=0 

where the n coefficients are functions of the system 
parameters (r=l,2, e . « ,m) . 

A variation of system parameters q will change the coefficients a 

r ic 

resulting in a change in the n root locations. In order to evaluate 

the change in the root locations, the real and complex roots are 

considered separately. 

th 

Consider the i pair of complex roots of equation (7-1) as: 



’^i.i+l 




(7-2) 
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til 

and the J real root of equation (7-1) as: 






(7-3) 



Macroscopic sensitivities, the sensitivities of the roots (7-2) and (7-3) 

V 

to large (finite) changes in parameter q^, are defined as: 



i,r 






q^-0,li^r 



l,r 



Ato 

j 



q^-0,l#r 



(7-4) 



J,r 



Ao 

Aq 



q^-0,l/r 



S. , macroscopic damping sensitivity, is a measure of the change In 
1 , r 

the damping ratio of the complex root pair r^ due to a finite change 

‘^i ■; 

in parameter q . S , macroscopic natural frequency sensitivity, is 

IT X I IT i 

a measure of the change in the undamped natural frequency lo of the 
complex root pair r^ due to a finite in parameter q^, Sj , macroscopic 
real root sensitivity. Is a measure of the change in the real root r^ 
due to a finite change in parameter q , 



■J 



Microscopic sensitivities, the sensitivities of the roots (7-2) 
and (7-j3) to small (infinitesimal) changes in parameter are defined 
as : 



8 



S 



l,r 

0) 

"i 

l,r 






3o) /3q 

n, ‘^r 



(7-5) 



s. - 9o./3q 

J.r J 
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The evaluation of microscopic sensitivities equations, (7-5) is 
described in detail in reference [14] and is reviewed in reference [5]. 

In essence it consists of writing the characteristic polynomial as two 
equations in the Chebyshev functions of the second kind, equations (2-6); 
computing the partial differentials the functions with respect to the 
parameter q^, and o)^; rearranging the results to get two simultaneous 

equations in and with coefficients and constants which 

are summations in and ^ and a^; and, solving for the 

desired sensitivities. 

Macroscopic Sensitivity of Singular Lines 

The sensitivity of singular lines to finite changes in parameters 
a and 3 is determined by applying equations (7-4) using values of a, 3, 

C, and 0 )^ obtained from the parameter plane diagram. 

Example VI 

Consider the characteristic equation of Example V and the para- 
meter plane diagram. Figure 4. Assume an operating point (25, 9.3) 
on singular line 23 for which c = 0*3 and = 12.02. To move to 
singular line 33 involves an infinity of complex root sensitivities 
since the singular complex roots are represented by an Infinity of points 
on line 33; however, the following are considered to be of primary 
interest : 

1. Movement of the operating point on a line normal to the terminal 
singular line. In this case movement to point M^ (24.6, 12.7) on 
line 33 for which ^ = 0.4 and = 12.37. The macroscopic root 
sensitivities are 
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0.35 

24.6-25 



0.1 

24.6-25 



-0.25 




a 



0.35 

-0.40 



-0.875 



0.1 

12.7-9.3 



b) 



= +0.032 



0.35 _ 0.35 

12.7-9.3 ■ 3.4 



= +0.103 



2. Parameter S fixed. In this case the operating point moves to point 
(17.7, 9.3) and the macroscopic root sensitivities are 



• = -0.013 S " = = -0.048 

a 17.7-25 a -7.3 



3. Parameter a fixed. The operating point moves to point (25, 12.9) 
and the macroscopic root sensitivities are 




a a 



S 



4 

6 



0.1 

12.9-9.3 



+0.028 



“n _ 0.35 
6 ■ 3.6 



0.097 



The sensitivities computed in Example VI show that the frequency 
of the singular line complex roots is more sensitive than the damping 
ratio to finite changes in both parameters singly or together. Inspec- 
tion of Figure 4 shows that this result is true for all singular lines 
in the stable area since they are uniformly divergent with increasing 
positive slope for larger damping ratio and smaller natural frequency. 
In addition. Figure 4 shows that the sensitivity of singular complex 
roots to parameter variations is decreased as the operating values of 
a, 8, or a and 8 are increased; that is, as the singular lines diverge 
sensitivity is decreased. 
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Complex root sensitivities of the singular lines could be used 
effectively in self adaptive control systems. For example, consider a 
system with characteristic equation (5-2) for which the desired operating 
point is M2(29. 6,12. 7) on singular line 33 of Figure 4. Consider that 
either parameter a or 3 changes from its design value due to aging, 
breakdown, or external disturbance resulting in changes of both C and 
0 )^ of the complex roots of the system. As shown in Example VI, natural 
frequency is more sensitive than the damping ratio to finite changes in 
both a and 3. A sensor at the output could measure the system response 
to a standard input, compare this with a model generated output, and 
use the difference in transient frequencies to determine the change in 
natural damping of the dominant roots. A controller could then determine 
and apply the change in parameter 3 necessary to drive the operating 
point back to singular line 33. Since the complex roots are dominant, 
the real root near the origin is fixed, and the remaining real root is 
far removed in the left half s-plane for any point on the singular line, 
the specific point of return to line 33 is not significant. Parameter 
3 is therefore chosen as the correction element since 

0 ) 03 

^^e^Aa=0 ^ ^^6 ^Aa=0 ^ ^A3=0 
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VIII 



APPLICATION OF SINGULAR LINE THEORY 
TO LINEAR CONTROL SYSTEMS 



In this section selected linear multivariable control systems are 
considered in sufficient detail to demonstrate the applicability and 
potential of singular line theory. The analysis of each system is not 
complete; discussion is limited primarily to the singular lines on the 
parameter plane diagram and their interpretation. 

The basic system to be analyzed is an inertially stabilized space 
vehicle. [13] Specifically, the system is comprised of two, cross 
coupled channels each with reference and disturbance inputs, plant, lead 
compensator, coupling element, and unity feedback path. Throughout the 
succeeding examples complete symmetry of plants, compensators, and 
coupling coefficients is assumed. That is: 




G « G = G 
Cl C 2 c 



( 8 - 1 ) 

^ 11 ° ^22 
^12“ ^21 

Examination of the system flow graph, Figure 5, yields the 
characteristic equation: 

A » l+2a,-G G+(a,,La 2)(G G)^ (8-2) 

11 c 11 Iz c 

Example VI 

Assume plant and compensator transfer functions: 
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Linear, Multivariable System Flowgraph 



G 



2 

s 



G 

c 



=- (■ 



TS+1 

YTS+l 



) 



(8-3) 



Substitution of equations (8-3) in equation (8-2) and definition of para- 

2 2 2 

meters q^ = 2a^^KK^ and q^ = (^2.1 ”^12 ^ gives the characteristic 

equation as: 



(yt)^s^ + 2yts^ + ('^+Y'^)Q2^s^ + 



2 2 

(qf+T q2)s + 2xq2S + q^ = 0 



(8-4) 



Using the approximate analytical methods of reference [13] 
compensator parameter values of t = 0.5 and y “ 0.05 are chosen to 
provide dominant roots with C = 0.5 and ^^=1* Consider parameters q^ 
and q^ as variables a and B in the parameter plane coefficient equation, 
equation (2-9). The characteristic equation then becomes 



A = s^ + 80s^ + (l600+20a)s^ + 840as^ + 



(1600a+4003)s + 1600Bs + 1600B = 0 

The system transfer functions are 



(8-5) 
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The numerators of equations ( 8 - 6 ), (8-7), and (8-9) are real functions 

of one relatively simple parameter formed from the gains and coupling 

2 

coefficients of the system, i.e., ^2/q^, ^12^1/a^^, and ^12 ^1/2. The 

zeros of these transfer functions can, therefore, be obtained by root 
locus methods. The numerator of equation ( 8 - 6 ) is a function of para- 
meters K and q^^; the zeros of this transfer function can be specified 
by two methods : 

1. Fix K and vary q^ to determine the root locus, 

2, Designate q = a, *7 = 3 and use the quadratic coefficients parameter 
plane method to obtain a parameter plane diagram. 

These transfer functions are included simply to indicate a possible 
analysis method; no discussion of them will be given since it is outside 
the scope of this study. 

Figure 6 is the complete parameter plane diagram of equation (8-5) 
for selected values of C and co^. Associated computer output data is 
contained in Appendix V, Table V. The constant C and constant curves 

in Figure 6 were predicted by previous parameter plane theory, but the 
singular lines were not. 

For system specifications requiring an underdamped or critically 
damped response, previous parameter plane theory restricted values of 
parameters a and 3 to the section of the parameter plane bounded by the 
a = 0 contour and the C = 1 contour. Thus values of a and 3 which 
specified an operating point in the stable real root area were prohibited. 
This example shows, however, that the singular lines define operating 
points at which complex roots exist throughout the area previously 
considered solely a real root area. Each singular line is tangent to an 
associated constant C curve. The value of co on the constant ? curve at 
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the point of tangency defines the constant ^ and constant cd of the 

n 

singular line. In addition, in this particular case, the parameter 
plane diagram shows that: 

1. Each singular line is not only a line of constant C - constant o) 

n 

for one pair of complex roots, but is also a line for which the 
value of one real root is constant. Thus the singular lines are 
coincident with constant a lines. Specification of the two system 
parameters and such that they define any point on a singular 
line fixes three of the six roots of equation (8-5). The remaining 
three roots vary with the position along the singular line. For 
example, singular line 101 for which ^ “ 0.9 and o)^ = 3.83 has 
fixed roots: 

®1 2 “ -3.448 - j 1.670 
= -33.103 

2. The intersection of two singular lines defines the values of the 
system variable parameters q^ and q^ necessary to fix all six roots 
of equation (8-5). For example, operating point M^(64.5, 650) which 
is at the intersection of singular lines 72 and 101 results in 
characteristic equation roots: 

Sl/2 “ -33.92 - J25.22 

= -21.648 determined by line 72 and, 

= -3.448 - jl.67 

s. = -33.103 determined by line 101. 

6 

This was to be expected since such an intersection defines unique 
values of a and 3 so that the coefficients of equation (8-5) are 
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all constants. Previous theory predicted the intersection of six 
constant a lines rather than two singular lines each coincident with 
a constant a line and, in fact, six intersecting constant a lines do 
not exist. 

3. Singular line 41 for which ^ = 0.3 and = 1.26 clearly defines a 
dominant root line in the stable area of the parameter plane. 

Despite variations in parameter a and 3 along the line, the roots 
s = -0.3793 - j 1.2062 remain dominant, the real root o = -39.24 
remains fixed, and the other roots move to satisfy the characteristic 
equation. Inspection of the computer data printout for singular line 
41, Table VI of Appendix V shows that the singular roots have a 
dominancy of about three at point M2(10,13.2), but as the operating 
point moves along the line to M^(80, 123) the dominancy is reduced 

to approximately one point six. 

4. Theoretically a dominant root area bounded by chosen singular lines 
diverging from singular line 31 could be defined. The choice of 
specific lines would be determined by the degree of dominance 
acceptable; for example, the singular roots have dominancy of between 
2.3 and 1.2 as singular line 31 is traversed. This example shows 
that the approximate analytical method of reference 13 is at best 
very approximate. The required dominant roots were C = 0.5 and 0)^=1 
with a dominancy of about ten, but, in fact, on singular line 51 with 
C = 0.5 and o)^ = 1.055 the dominancy is between 1.5 and 1. 

Figure 7, Sensitivity of Singular Lines to Parameter Variations, 
is the parameter plane diagram of equation (8-5) for values of C between 
0.5 and 0.6 and C between 0.8 and 0.9. Inspection of the diagram yields 
the following information. Maximum sensitivities of C and of the 



57 



o\o 0ie on oiG oi^ oea occ 




K iCfLl • J.e8e^8) UNJT3/Jti,H 

'i-jcite > o.6el£^e^ UKITS^JMCH. FIGURE 6 

SINGULAR LINES ON THE PARAMETER PLANE 
OF A MULTIUARIABLE , COUPLED SYSTEM 



58 



singular complex roots to finite variations in a, 3, and a and 3 occur 
at the points of tangency of the singular lines with associated constant 
C contours. That is, points (8.6,22) for C between 0.5 and 0.6 and 
M 2 (2.175, 119) for C between 0.8 and 0.9. As singular lines diverge 
from these points the macroscopic singular complex root sensitivities, 
equations (7-4), decrease. 

Consider operating point (36.35, 135.6) on singular line 11 

for which c = 0.5 and o)^ = 2.1115. Moving the operating point to 

M. (36.35, 141.2) on singular line 21, c = 0.51, o) = 2.154 yields 
4 n 

sensitivities : 



Aa=0 



0.51-0.50 _ 0.01 

141.2-135.6 5.6 



0.0018 



0) 




Aa=0 



2.154-2.1115 _ 0.039 
5.6 ^ 5.6 



0.0070 



That is, the natural frequency is approximately four times more sensitive 
to changes in 3 with a constant than the damping ratio. Moving the 
operating point to M^(35.3,135.6) on singular line 21 yields sensitivities 







0.01 


0.01 


a ) 


A6=0 


" 35.3-36.35 


-1.05 




1 


0.039 




“ J 


1 

A 6=0 


-1.05 





-0,0095 



-0.037 



That is, the natural frequency is approximately four times more sensitive 
than the damping ratio to changes in a with 3 constant. Moreover, a 
finite change in a has approximately five times the effect of an equal 
change in 3 on the natural frequency and damping ratio of the singular 
lines , 
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Example VIII 



Consider that the specification for the basic system’s response 
is changed from underdamped with one overshoot to critical damping. 
Assume the same system structure and transfer functions. The system 
characteristic equation is equation (8-4). The approximate analytical 
methods of reference 13 now require compensator parameters t = 0.707 and 
Y = 0.05 to place dominant roots at C = 0.707 and = 1. Considering 
parameters q^^ and q 2 defined in Example VII as variables a and 3, the 
system characteristic equation is 



The parameter plane diagram for equation (8-10) is shown in 
Figure 8. Computed data for the singular lines is contained in Appendix 
V, Table VII. The constant C curves are similar to those in Figure 6 
or Figures 2 and 3 in that they divide the stable region of the parameter 
plane into two distinct areas. The stable region is bounded by the C = 0 
contour and the positive a axis. The section of the stable region 
containing the intersecting constant C curves was previously considered 
the complex root region; existing parameter plane methods correctly 
solved for all six roots of equation (8-10) in this region. The section 
of the stable region bounded by the C = 1 contour and the positive a axis 
was considered the stable real root area. Singular lines exist in this 
area. Previous parameter plane methods determined either four or two 
real roots for each point in the stable real root area, displaying the 
reaults as intersections of constant a lines. Inclusion of singular 



! 



+ 56.56s^ + (20a+800)s^ + 593. 9as^ + 
(800a+4003)s^ + 1131.56s + 8006 



0 



(8-10) 
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lines completes the parameter plane diagram allowing graphical solution 

for all roots of equation ( 8 - 10 ). 

Figure 8 , the parameter plane diagram, shows that: 

1. The value of constant for a singular line is determined by the 
value of on the constant c curve at the point of tangency of the 
singular line. For example, singular line 51 (C*0,5, a)^*1.4929) is 
tangent to the 5-0,5 curve at point Mj^(3.95, 3.91) whereas line 52 
is tangent to the 5 - 0.5 curve at point M2(76.4, 1460.). 

2. Each singular line is coincident with a constant o line. For 

example, line 71 for which 5-0.7 and o) -2.0972 has fixed roots: 

n 

8 ^ 2 “ -1'^68 - J 1.4977 
83 - -25.34 

3. As in Example Vll, the intersection of two singular lines define the 
value of parameters a and 8 necessary to fix all roots of equation 
( 8 - 10 ). 

4. The singular complex roots become dominant at singular line 61, that 
is for 5*0.6. The amount of domlnancy is Increased as the operating 
point moves to singular lines 71, 81, and 91, i.e., as 5 Is 
increased to 0.9, but the maximum achieved is approximately 1,9. 

5. The natural frequencies and damping ratios of singular lines 11 
through 101 are more sensitive to finite changes In parameter 8 

than to changes in a. Since these singular lines diverge, effectively 
from the origin of Figure 8 , the singular root sensitivities increase 
with respect to a but decrease with respect to 8 as the operating 
point 8 value Increases with a held constant. 
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Example IX 



Consider the basic coupled system structure of Figure 5 but 
assume plant transfer functions: 

G - ■ ■ V ^ ■ (8-11) 

(s +2S+2) 

Assume the compensator transfer functions and parameter values of 
Example VII. In addition, consider system parameters and q 2 as 
defined In Example VII as parameter plane variables a and The 
resulting characteristic equation is 

s^ + 84s^ + (20a+1928)s^ + (880o+7048)s^ + (3320o+400@+13444)s^ + 

(4880a+16003+13120)s + (3200a+1600g+1600) - 0 (8-12) 

Figure 9 Is the parameter plane diagram of equation (7-12) for 
selected values of ^ and U)^. Table VIII of Appendix V, contains the 
computed data for the singular lines. As o)^ varies from 0 to » Che 
constant 4 contours curl around Che origin through the second, third and 
fourth quadrant of the a3 plane, reverse their direction producing 
sharp spikes In the lower section of the first quadrant, and then proceed 
to infinity via the first quadrant of the oi3 plane. For clarity only the 
final sections of these contours have been shown In Figure 9; the 
sections not shown are similar to the spikes of the 4*1 contour. The 
4-1 contour Is a particularly confusing contour even with small scale 
graphs and complete data printouts available. For values of u)^ from 0.2 
to 2.06 It describes short arcs in all four quadrants of the a3 plane 
with Intervening discontinuous jumps. For values of o)^ greater than 
2.06 the curve Is as shown In Figure 9 crossing the origin and proceeding 
to Infinity via the second quadrant. 
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The bounds of the stable root area are the C = 0 contour and a 



line determined by setting s=0 in equation (8-12), that is: 

3200a + B + 1600 = 0 

B = -(2a + 1) (8-13) 

Within this stable area there are four distinct sections corresponding 
to different real and complex root combinations. Previous parameter 
plane methods solved for three pairs of complex roots, two real roots, 
four real roots, or two real roots and a complex pair in these sections 
of the stable area. The inclusion of singular lines in the parameter 
plane diagram defines the remaining roots and allows graphical solution 
for all six roots of equation (8-12) at any point in the stable area. 

The parameter plane diagram and associated singular line data 
show that : 

1. Singular line values are accurately determined for specified values 

of Z and 0 )^. For example, line 51 for which ^ = 0.6, = 33.272; 

the complex roots s = -19.963 - j 26.618 are correct to three 
decimal places. 

2. No singular line can be considered as a dominant root line and 
therefore a dominant root area cannot be defined. 

3. The singular lines are technically not coincident with constant a 
lines. One spike of the C = 1 contour for from 1.93 to 2.26 
produces a mesh of double-valued, constant a lines which encompasses 
the stable real root area. The result is that, effectively, each 
singular line is coincident with a double-valued constant o line. 

4. In part of the stable area singular line intersections occur solving 
for four roots which were previously undetermined. 
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5. When a singular line is extended into the fourth quadrant of the a6 
plane below the stable area boundary defined by equation (8-13) the 
complex roots defined by the singular line remain fixed but the 
other roots shift to produce a right half s-plane root. For 
example, since line 73 crosses the stable area boundary as a is 
decreased from 40.0 to 20.0 a system root moves from o = -2.02 to 
0 = +4.97. 

Example X 

Consider the system structure, parameter values, and parameter 
plane variables of Example VI but now assume plant transfer functions: 



G - K 

(s+1) 

The resulting characteristic equation is 

s^ + (20a+42)s^ + (20a+4008+161)s^ + (2440a+16008+3280)s + 
(1600a+16008+1600) - 0 



(8-14) 



(8-15) 



The parameter plane diagram for equation (8-14), Figure 10, shows 
that no singular lines exist. Values of 4 between 0 and 1 were investi- 
gated. For all values of 4 less than 1 the coefficient matrix of the 
parameter plane solution equations, equation (4-4), is non-singular. 

For 4*1 there are two values of o) for which the coefficient matrix is 

n 

singular, but for both values the parameter plane solution equations 
(4-3) are inconsistent. 
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IX 



CONCLUSIONS AND RECOMMENDATIONS 
FOR FURTHER WORK 



Singular line theory is a small but meaningful addition to the 
parameter plane method. Singular lines, constant damping ratio - 
constant undamped natural frequency lines, occur in only a limited 
number of cases, but where applicable their addition to the parameter 
plane diagram defines previously undetermined roots. As a result, 
complete solution of all characteristic polynomials in which two 
variable parameters appear linearly in the coefficients is now possible. 

The conditions for existence of singular lines are complex. 

General rules for predicting singular line existence were developed, 
but ^n existence theorem which specifies the exact combination of 
coefficients that allows singular lines should, if possible, be developed. 
Such an existence theorem would facilitate synthesis and design of 
singular line systems. 

The computer program developed during the course of this study, 
although relatively slow, provides an accurate means of solving for and 
graphically displaying singular lines in the parameter plane diagram in 
conjunction with constant cu^, and a curves. More efficient program- 
ming techniques and the use of a faster polynomial root solving sub- 
routine should be considered. 

Dominance of singular roots results in a dominant root line on 
the parameter plane diagram. This allows choice of values of parameters 
a and 3 from the infinity of points specified by a singular line while 
simultaneously ensuring dominant mode operation of the system. 

Sensitivities of the damping ratio and undamped natural frequency 
of singular roots to finite changes in parameters a and 3 were computed 
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for specific linear, multivariable control systems. It is considered 
that measuring the movement of singular roots resulting from finite 
parameter variations could prove to be a useful concept in self-adaptive 
control systems. Since movement along a dominant root singular line 
implies little change in operating conditions, the adaptive controller 
need only drive the system back to the singular line not to a specific 
operating point. 

Singular line theory is based on the special case of linearly 
dependent parameter plane solution equations. The theory may not, 
therefore, be extended to include polynomials whose coefficients are 
linear functions of two parameters and their product. 
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APPENDIX I 



TABLE OF THE CHEBYSHEV FUNCTIONS OF THE SECOND KIND 



The Chebyshev functions of the second kind, are defined 

by the recursion relation: 






where: 



Uq(C) = 0 



u^(?) = 1 



These functions for k = -1, 10 are 



k 

-1 

0 

1 

2 

3 

A 

5 

6 

7 

8 
9 

10 



-1 

0 

1 

2 

- 1 

- Ac 

16c - 12c^ + 1 
32c^ - 32c^ + 6c 
6Ac^ - 80c^ + 2Ac^ - 1 
128c^ - 192c^ + 80c^ - 8c 
256c® - AA8c® + 2A0c^ - AOc^ + 1 
512c^ - 102Ac^ + 672c^ - 160c® + lOc 
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APPENDIX III 



EXPANSION OP 



S ^2 



Given a characteristic polsmomiil of the form; 



n . 

f(s) * ^ (b^o<+ cjjf + d^) s' 



k«o 



(A-1) 



Develop a general expression for: 



®2 ^2 



• ®1°2 ’ Vi 



(A-2) 



where ; 



2: (-o\ 



"2 ' 



ic n 

k*o 



(k) 



(A-3) 



!, .^(-1 )\ , Cj (-1 )\ 



k*o 



(k) 



k*o 



Consider n 2: 



Equation (A— 1 ) becomes; 



(b2(<+c2p «»-d2)s + (b^o^+c^p+d^s + (b^o(<.c^+d^) , q 



74 



Evaluating equations (A-3) yields : 



B, = 



®2 ■ -‘’I "n + '■ 2 "^ 



C, . -c^ + Oj 



C_ -c. O + c„ 6) (2 5 ) 

2 In 2 n 



From which: 



'^iWl ‘ Vl"n-Vl“n^-V2"n^<2^> - * 



+ ‘'l'=2"n^ “ ° 



After grouping terms, this is : 



B,C2-B2C, = 



(boC,-b,Co)«„ - (V2-Vo)"n 



(A-4) 



Consider n 3 : 

Equation (A-1 ) gives : 

(b^o^+Cj p +d^)s^ + (b^*^ +c^p +d2>s^ + (b^ ^ ’'■‘^1 '*’‘^ 0 ^ * ^ 

Evaluating equations (^ 1 .- 3 ) yields : 

j 3 = -b + b^«o ^ - b_co ^(2^) 

1 o 2 n 5 n ' 

\ "n ^ ^2"n^^2^ ) - b^w^5(4S^-1) 

S • ■‘^0 •" 

C2 • -c^ co/(2 S ) - C3 ^/(4 X^-1 ) 
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iMOte that in each case only one term of the summation need be calculated, 



i.e; 



1 



(n) 



d , + (-1 )\ 1 ^ 

V-1) "" “ 



(A-5) 



Mviltiplying as indicated in equation (A-5)» cancelling and grouping 
terms gives : 

“ ^V2-^2°o^ ) + (b^c^-b^c^)^)J(4S^-l) 

a, +(b^C2-b2C^)«J - (b^c^-b^c^)«^'^(2t ) 0 



Mote that only the tenris involving b^ and c^ are new, that is : 



(B, 






(n) 



(n-1) i»o 



Consider n 4 : 



evaluating equations (A-1 ) in accordance with (A-5) yields 



B » B 

^ \nO) ^ " 

B » B +b o/(8X^-4X) 

^ ^(n.3) ^ " 

C r C •►C «AsX^-45) 

^ ^(n»3) ^ '' 
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Substituting in equation (A-2), cancelling and grouping terms gives : 



+ • (k.Cj- t (iX'l) 

+ - (kS"^40^Vax) 



Extension to the rreneral case ; 



Vx-BxC, . 






•+ [^f.i'v-i»jCf.,](-ir4-M<' * 







“t c 

U A-l /K 



a^- 1 



where 


i 


1,2,... 




is 


the 


roYT nixmber 




J 


1,2,... 




is 


the 


column number 



(A-6) 
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Or, in compact notation ; 



“IW) 


= z 

i.j 


where : 


CE * 


L ij] 



ce. , 



(W) 



whose elements are : 




c 



i-f-j-1 



- b. . , 







2i+j-2 

n 



i+j < n-2 
itj ^ n-2 



(A-8; 
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APPENDIX IV 



PROGRAM PARAM S 



••JOB0629F» BOWIE D J 000 

PROGRAM PARAH S 001 

C THIS PROGRAM IS APPLICABLE TO POLYNOMIALS WHOSE COEFFICIENTS ARE OF THE 
C F(M*M (B* ALPHA + C*BETA + D) WHERE ALPHA AND BETA ARE VARIABLE PARAMETERS 

C AND 6» C# AND D ARE CONSTANTS. THIRD PARAMETERS CAN ALSO BE SPECIFIED 
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PRINT 220 140 

527 FORMAT (//*15X*5HALPHA*16X*4HBETA*15X*5H0MEGA*16X»4H2ETA*5X» 1410 

115HTHIRD PARAMETER.//) 142 
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702 IF( IPRINT)459.459.223 250 
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TABLE IV 

ROOTS OF THE CHARACTERISTIC POLYNOMIAL AT SELECTED POINTS 
ON SINGULAR LINE 43 OF FIGURE 4 



UJ 

z 

QC 

< 

Ol 



o 

o 



o 

o 

o 

o 



^fSiCsJO orsJCNjO' rocsjcNjo Ofsjrvio 

^ooo oooo ^ooo oooo oooo 

LLILLIUJUJ LULiJUJ UJ UJUJUJUJ UJUJUJUJ LU LliLUUJ 

tofvjrgo >oooo oromo 

H-r^cgfVJOH-ooa"OKOsjv/>v/>OH-ocDaoo^-ooa'^o 
qc pH »-i o QC o ^ o in a: eg O' o oc o O' O' o QC o O' cMn 
<ir\r-4.-^o<ooo>o<^ooo<oooo<ooo>o 

II 





o 


>• 


1 1 > 


1 >- 


1 1 > 


1 1 > 


o 


• 


QC 


QC 


or 


VC 


VC 


oc 




< 


< 


< 


< 


< 






z 


z 


z 


z 




I 






k-H 












<3 


o 


o 


o 


o 






< 


< 


<x 


< 


< 






X 


X 


X 


X 


X 





o 


00 


OrHrHOg 


Of-^.-4og 


o^«^ro 


Oi-4^fO 


3^i-4r0 








OOOO 


OOOO 


oooo 


oooo 


oooo 


UJ 


LU 


a 


g* ♦ 






4- 4> 


4- 4- 4- 4 


a 


in 


o 


UJUJUJUJ 


UJLUlULU 


UJUJUJUJ 


UJUJUJUJ 


UJUJUJUJ 


o 


fO 


oc 


00 <N eg in 


eginmoo 


ro 00 00 in 


eg4‘^'0 


oo»^r^r^ 




>0 




•^egegro 


^r-r-o 


rgmen^ 







00 



a 

UJ 

o 

C£ 

UJ 



CM 

>C 



O 



h- o GO 00 m H* r* CT' O' ^ rr> %o vO 00 >o o CO h- o o 00 
0^ r- fo m ^ 0^ 0"<r ^ 00 0^ eg m ro CNJ QC ^ fo ro rg 0^ ^ rom fvj 
<e- 0'Oe-<oog3oeg<roooNO^<oo>o>oeg< co^oom 

• • • ♦ 

I I I I 



CL • • 


• «cv 


• • ♦ •Ql • • 


• •CL • 


• • 


1 1 


1 


1 1 1 1 II 


1 1 1 


1 1 1 


«J 




•J 






< 


< 


< 


< 


< 


UJ 


UJ 


UJ 


LU 


UJ 


(C 


Qc: 


oc 


OC 


oc 



2T 

k-4 00 

in 



•-I • 



X 


1 






o 


*H 




eg 


< 






o 


o 


o 


o 


O 


X 






4- 


4 


4 


4 


4 






< 


LU 


LU 


UJ 


UJ 


LU 






H- 


ro 


O 


eg 


o 


O 






LU 


h- 




o 


oo 


O' 






00 


m 




eg 


in 






o 




pH 


CT' 


#H 


eg 


o 




o 




eg 


cr 






pH 




4- 




• 


• 


• 


• 


• 



UJ I 

< o 

H- O 

LU O 

K4 O 

in 



< 

o 

LU 

z 

o 










99 



20000E+02 .13521E+02 REAL PART IMAGINARY 



(\J<MO OCSJCVJO OCNJCsiO OCSjrVjO 

F^OOO oooo oooo oooo 

LULLMJUUJ UJUJUJLU LULUUiULI UiUiUJUJ 

m^O'O ornroo OtfMAo or^r^o 
h- fO 0^0 OH- oooo K-OOOO)- oooo 
arHO^o^OQ^ooooQCooooaoooo 
<<MOOO<0<-*t^O<OHiMO<OH^O 
Q, ^ o Qu o ^ O G. o «-i H O O. O H ^ O 

#••• #••• 

I f!> lt> lf> It 
Ct u CC 

< < < 

z z z 



o 

< 



o 

< 



o 



OHFHro o f-4#hco o> <fO 

oooo oooo oooo oooo 

-f ♦ ♦ ♦ 

UJUJUJUJ LULULiJUJ LLIUJUiUJ UJUJUJUJ 
F-iOOh* cgCNOJH 

O ro fo f\i a o rn CO <\j oc O ro fO fNj QC O rri rn f\i 
00 o < 00 ^ ^ m< 00 >c ^ o < 00 ^ H 

• • t *0. • • • #0. f • t tO. • • f • 





«J 




< 


< 


< 


UJ 


UJ 


UJ 


cc 


cc 


a 



eg 


eg 


(M 


O 


O 


O 


4- 


•f 


4 


UJ 


UJ 


UJ 


fO 


ur\ 


H 




00 


•M 


>0 


H 


O' 




O' 


INI 


**4 


Fg 


INI 


• 


t 


• 



<NJ 


eg 


eg 


O 


o 


o 


♦ 


♦ 


♦ 


LU 


UJ 


UJ 


O 


o 


o 


o 


o 


o 


o 


o 


o 


m 


o 


in 


eg 


fo 


rn 


• 


• 


• 
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COMPUTED DATA FOR SINGULAR LINES OF FIGURE 6 



UJ oooo oooo OOOOOOOOOO o 

oooo oooo OOOOOOOOOOO 

3 u ULiULIUJULi UUUJUUUJ UJUJUJUJUUIJUIJUIJLIUJIJLIUJ 

<l OOOO OOOO OOOOOOOOOOO 

OC OOOO OOOO OOOOOOOOOOO 

< oooo oooo OOOOOOOOOOO 

Q. OOOO OOOO OOOOOOOOOOO 

OOOO OOOO OOOOOOOOOO o 

o •••••••••••••••••«• 

QC 



o 

o 

•f 

UJ 

o 

o 

o 

o 

o 



UJ 

a 

o 

•t 

c/> 



o 

UJ 



00 

K 

< 

X 



< 

H- 

UJ 

fSI 



< 

o 

UJ 

X. 

a 



oooo ooooooooooooooo 

UJUJUJUJ UJUJUJUJ UJUJuJUJUJLUUJUJUJUJUJ 
ir> 4 'oaom<\j< 7 'f\i o'o*-<rooc 7 'a>o cj> 
O'QOO'O Gooom fvjh-r^o^ocra''^m 4 » 
fM^-' 4 *co®' 0 ’ 00 *-<rOro 
r^ao oi/> ioor^*-H rsjincT'roooa^fMfNjo^’ 



0 ' 4 ' 0 rn*-<m 

OO OOOOOO OOOOOOOOOOO 
UJUJ UJUJ UJLU UJUJ UJUJ UJUJUJUJ UJUJUJUJ UJ 

ooomcsio^fvj c^^-c^Or-i^noc^a'c^c^ 
O'GocDOoooom <>ar^r»-o-<oa'a'' 4 “ro>r 
ifSsO >tQ 0 G 0 ^% 0 O'-Hr 0 m 

oirmo fvi goo c^ eg rg^ 4 - 

rg*-< irkinr^' 4 ’C 7 ^f 0 ^rg *-4 



00 OOOOOO OOOOOOOOOO-^ 
00 OOOOOO OOOOOOOOOOO 

UJUJ UJUJ UJUJ UJUJ UJ UJUJUJUJ UJUJ UJUJ UJUJ 
00 OOOOOO OOOOOOOOOOO 

0000000000000000000 

00 OOOOOO OOOOOOOOOOO 
00 OOOOOO OOOOOOOOOOO 

rgegrom mmosor^r*-aoooa'0'— < 



ornofNJ— «fNJ*«irg 

00000000 OOOOOOOOOOO 

UJUJUJUJ 111 UJ UJUJ UJUJ UJUJ UJUJ UJUJ UJUJ UJ 
oooGOGo^o^fnm 

'j'ln 4” nO r^— < — ^GOf^cginr^som— <G0-^ 
rHO' cmO' ocg oo'4‘ *hoo fom^o^o^mroGOh- 
ogGO rgroor^ ^r^in—^O'^omrOGdocM 
nJ-fh caocT' ^> 4 ' (Nfn f\jfnf\if\jrnrgfOf\j^ 





. 22301 E +02 


• • 

z z 
0 0 


M »— t 


eg 


H- h* 


0 


0 0 


4 - 


«J 


UJ 


Q 0 




00 00 


0 

rr\ 


0 0 


eg 


Z Z 


eg 

• 


# • 

z z 

UJ UJ 
h- H 

00 00 






OOOOOO 


0 


+ z •♦■z 


4 - 


UJOUJO 


UJ 


oooo 


0 


ozoz 


0 


0 ►— i 0 


0 


0 0 


0 


•- 4 O 0 — <oo 


r -4 


•z «z 


• 


0 0 

H H- 
< <l 

0 0 

0 CJf 
UJ UJ 




eg 


OLUOIL 


0 


-f-O -♦-o 




UJ UJ 


UJ 


ozoz 


GO 


oujouj 


eg 


Oh- Oh- 




000 000 


GO 


cg>-rg> 




•00 *oo 


• 



•-irg— —<rgpHrg^— <cg*- 4 rg#-^<\j.^c>g^rg#-i 4 
rgrgmm ^^j-inin 'OOr^^-oocoo^a^oo— < 

fP^#P^^n4 

^gp^^g^g kjnj Mfsj ^g^g ^g^gKl^sl^gKJ^g^gNJ 
ISIS ^:si:s 
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TABLE VI 

ROOTS OF THE CHARACTERISTIC POLYNOMIAL AT SELECTED POINTS 
ON SINGULAR LINE 41 OF FIGURE 6 



LU f-H 

Ol O 

O 4- 

UJ 

</> 4TN 

00 
>o 
in 



oooooo 
I I + I I 
UULPLUUJUJUJ 



•■-4 p-H rM f— < c 7 ' cr* 
oooooo 

^ + I I 

UJUJOJOJUJUJ 



OOrHOOO 
+ I I I I 

LUOJUJUJLUUJ 



O' (\^ J 
rHOOOOO 

UJUJUJUJUJUJ 



csjcmo'CT'O*-* rvj<\ioo>t(Njm >r fNifvcT'COoo 

fsj (\j f«i O' ^ >o mm >t O' h- vO sO >t o 00 ro »— m >0 >o >o O' 
oc ^ c^ o >c rn q: o o rH rH >r o oc o o m 00 nO cr >o o o ^ ^ 

<c m >o cNJ f\j rg nO < c\j (\j f-« ^ r- O' < f\j CM o nO nC < ro (Nj cNJ .-4 *>-4 



Cl--< 




•r? 0. •-"< 


r»* r— 4 rTj (\J CL rO »~4 1 — i r-H 


> 


\ \ 1 > 1 1 


1 > 


II 1 > II 1 


o: 


cr 


cr 


cr 


<T 


< 


< 


< 


2T 


2L 








l-H 






o 


o 


o 


O 


<3 


< 


< 


<l 


2. 


2- 


2. 


2. 


M 






1— 1 



CL 

UJ 

o 

(T 

UJ 



C/) 



o 

+ 

UJ 

m 

00 

nO 

m 



‘✓5 

»— 

U 

u 



iO':^r\irv 

OOOOOO 

UJUJLULULLJUj 
corvjcNjCT' rsj»~t 
o^o»-4mro^ 
h~r*vOroinr-rvj 
CT(M>Or-COOC7^ 



oOrHf-irvjoj oo-H.~4rsi(\i 

OOOOOO oooooo 

UJ UJ UJ UJ UJ U.J UJUJUJUJUJUJ 

(y>irr|,- 4 ,-^r -40 rO r<> O 

rO(^oo>t<\J fomoo^(\i'4’ 

K- cr cr ^r <\j m H- ( 7 ^ cr 00 O' vO fvj »— or' cr (\i 00 00 

cr r- fsj rvj cr m cr r-- r- r- m r- cr cr fTi r- r- cr oc 00 



OfMCVi^NJ 

oooooo 

•f -f -f -f -r -f 

UJUJUJUJUJUJ 
m ro rr> rH 'j- nJ- 
(\JcrifO>3‘^ 



< cr ^ rn st <r<'fr (\j(\jrO(^<l r<^ r^f\iCr<\;c^<rvfr 



CL • • • 

1 1 


III III 


III ill 


1 1 1 


1 1 


III III 

-j 


III III 


1 1 1 

— 1 


< 


c 


<T 


< 


LU 


LU 


LU 


o« 


cr 


cr 


cr 


cr 



X 

< 

X 



< 

K 

UJ 

M 



O 

o 

+ 

UJ 

o 

o 

o 

o 

rO 





^4 


(M 


(M 


CM 




O 


r, 


O 


O 




-f 


+ 


+ 


•f 


c 


LU 


UJ 


lu 


o. 




rg 


u^ 


o 


m 


IL 


C 


CM 


r— 4 




X- 


nO 


r 


cr 


m 




nT 




oc 


>r 




CM 


•—4 


CM 


>r 




• 


• 


• 


• 



I 









O 


(M 


CM 


(M 














4. J 


<L 




< 


+ 


+ 


+ 


+ 


O 


o 


X 


\±i 


X 


X 


LU 


Ux 


-f 


CL 


C]j 


o 


o 


o 




LU 




o 


o 


o 


Ct 


o 




<3 


o 


o 


o 


o 




'4’ 




o 


o 




o 




>o 




'O 




CM 


m 




CM 




• 


• 


• 


• 



fsl 

je 
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4COOUE+J2 .6028CE+J? REAL PART IMAGINARY PART 



C'gC'g 

oooooo 

UJUJU.< LULU'JJ 
'^rn(NQOvOvn 
mso>oooo(7'>— 
— «ooo'o<-)a: 

OC\iCsiNOOOCD< 

ij »«i-J ^"\J ^—4 ^ 



oooooo 

I + -f I «f -f 

LL’LiJUJiJJaJUJ 

I— 

OnOnOCOQOOO^ 

ocoocnjo^o^q: 

(NJCNJCvist fN(NJ< 



OOOOOO 

UJLUUJUJULuJ 
LO rn f\J fv 



OOOOOO 
I -f -K I -f + 

UJOJUJUJUJUJ 

om,-4^coQO 



oooooo 

I -f ^ I 4* 

LUVUUJUJULJIJU 



\Q sQ \Q(^ O Oh»rn' 0 '> 0 ' 0 Q 000 ^*“Q 0 ' 0'0 rsi *■ 
m o ooo o o q: rn o o <30 >r cx oo o (\i sO sO 
cNj (\jin r-. < rn c\j rsj in o O < rn f\l f\j in rn rn 

Csj Cu— <*—»^rnfnrn rrrr, 



> 1 1 


1 >- 1 1 


1 > 1 1 


\ >• 


Oi 


cx 


rX 


CX 


<r 


< 


<r 


< 


2T 


2: 




2L 




t— 4 


1—4 




O 


o 


o 


C5 


< 


<I 


C 


<1 


2: 






2. 


•—4 


fc-H 


•—4 





«-<OOrg<Njr..' 

OOOOOO 

LU LU LU UJ LU UJ 
rn r;~> I 
ao>T(Nj>th-r^ 

vTO^CT'njOCOO 

Ovjn-f^CT^cocox: 



^ rsj f\j 

OOOOOO 

-f 4- -f 4- -f -f 

UJUJUJLUUJUJ 
»— 4in 

v0>rCN4>i'OO 
00(T'0^(NJCrO>H- 
r-ih-fs-O' oocxDo: 



OOOOOO 
4-4-^-4- + 4- 

UJLULiJUJtJUlXJ 
r^rNiin.— iiT.ir. 
o>>r (\j>rro(N 
>r crcT'CNjO'cr 

1—4 r** r***CT' 00 00 



OOOOOO 
+ + + 4'4-4- 

UJLUUJUJUJUi 
O tXl ^ r-H 33 CO 
in>i'f\j^mrn 
♦— Osjcr O' fvicT' O' 
q^^^-n-o>oooo 



— 

OOOOOO 
4*4--f + 4«4* 

LULU LULU LULU 
mrvjun^sO'O 
h- <M >r nT 
i~OCT>o>rga>a' 
cx^f^r^O^oooo 



f\irnrnf<'^^<insjrnrnrn> i-4*-H<5(\irnrnrnp-i rnm»~4f-H<fvmrn rn#— 

• • • • • «a. • • • • • «CL • • • • • *0. • • • • • tQ. •••••• 

I I I I I I I I I I I i I I I I I I I I I I I I I I I I I I 
— I — I -j -j 



<r 

LU 

cx 



< <x < 

LU LU LU 

oc a: nc 



CN. 




m 


m 


r~'^ 


o 


o 


o 


4- 


-f 


4- 






lU 


LLl 


LU 


m 


o 


m 


nj 


nO 


m 


m 


O 


I.T 


vL- 


r- 


.-n 


tr^ 


r-H 


o 


(NJ 


r- 








• 


• 


• 


• 



(\i 


CxJ 


C\i 


CNJ 


I** 


o 




o 




4- 


4 


4 


LU 


LL 


LU 


LU 




C'\ 




O 


o 


O 


O 


CO 




O 


o 


o 




o 


« -J 


o 






h- 


X 


• 


• 


• 


• 
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TABLE VII 

COMPUTED DATA FOR SINGULAR LINES OF FIGURE 8 



or 



OI 


oo 


oo 


H- 


oo 


oo 


UJ 


+ -f 


+ 


X 


UJUU 


OJUJ 


< 


oo 


oo 


QC 


oo 


oo 


< 


C' * V 


oo 


Cl 


oo 


oo 




oo 


oo 


a 


• • 


• • 


oc 






1*^ 






X 















.-Hcn 


ocn 




oo 


oo 




1 + 


+ ^ 


Hi 


OJUJ 


UJUJ 


a. 


*— 4 in 


>or^ 


o 


-HO^ 


vCfsi 




> 4 - CO 


f-»r^ 




coo 


incM 




GOO' 


CO CM 




• • 


• • 



h- 






Q- 






UJ 


^CO 


Or<0 


O 


oo 


oo 


or 


1 + 


4 + 


LU 


LULU 


LULU 


H- 


OCO 


^vC 


21 


f-»CT' 


>OCM 




sTCO 






COO 


in CM 


00 


coa^ 


CO CM 




• • 


• • 


X 






< 






X 









oo 




OO 




oo 




oo 




+ + 




+ + 




LULU 




UJLU 


< 


OO 




oo 




oo 


• 


oo 


LU 


C30 


CM 


oo 


M 


OO 


• 


oo 






z-i- 


CM CM 




• • 


<z> 


• • 



xol xa 

Kl— V-h- 

uoo oou 

00 



LUUJ lULU 

«JCD -JCQ 





OC »1 




OCM 






oo 


LUK 


OO 


UJK 


< 


•f + 


CDO 




c:;u 


o 


LUUJ 


z 


LULU 


2 


UJ 


OCM 


K 


cn>o 


H* 


X 


>om 


O-J 


iTtcn 


O-J 


o 




Z-J 


in *-4 


20 




OCO 


f-H 


O'H- 


►-H 




rvji— » 


♦- JC 


U)vO 


H--S 




• • 


00 


• • 


00 






OH- 




OH 






xo 




xo 










o 




r-»CM 


ooa 


•-4 CM 


oOO. 






h- 




H- 




r ~4 


0 . 4 /) 


CM CM 


0.00 










x*-* 




*VJKJ 


OX 


MM 


ox 






sc:h- 


JK 3 E 


2 H 



oooooooooooooooo 

oooooooooooooooo 

UJ OJ UJUJ LIJ UJ UJLU OJ UJ UJUJ lU UJ OJUJ 

oooooooooooooooo 

oooooooooooooooo 

oooooooooooooooo 

oooooooooooooooo 

oooooooooooooooo 



o m f-» f\i *-• CM ^ (\i f\j 

CDOOOOO ->0000000 30 
+ + + + + + + + + + + + 
UJUJ UJUJ ULi UJ OIUJ OJ UJ OI UJ OJ UJ OJ UJ 
in 00 nO to CO o >o o r- nO ITS CX3 

CO ^ CO 00 ro fsj 00 o r- ^ r- »-• 

CO o 00 cr >c cr o O' >o o m 
r-- ^ m (\i CO Osi fsj m <N h- 



O *^nj fVJ r— tfSJ *— < (\J f\) •— < f\j r-*< (\j 

oooooooooooooooo 

+ + -*- + ‘f-*- + -f + + + + + + + 4- 

UU UJ UJ uU UJ UJ UJ UJ UJ UJ LLi tU tiJ UJ UJ LU 

^ o o rn r- O" >c o in f\j o (V 00 r- m 

v00>i-^mr->i‘a'^0'0'0r^cocor^ 
CO »-» c- ^ CO CO '4“ m CO o 
aoOrooo^r^O^NOC^ ooooro^^ 

.-I i-i uO C\l CO Csl f\l CO fvj *-» >0 0- *-• 



OOOOOOOOOOOOOO'-^*-^ 
OOOOOOOOOOOOOOOO 
+ + -♦- + + + + + + + + + + + 4- + 
IXJ IIJ LU UJ UJ UJ UJ UJ UJ LU LU UJ LU UJ UU LU 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 
OOOOOOOOOOOOOOOO 

m m ^ tn no >o 00 00 C 7 ^ O' 



o <VJ rvi •-» C\J C\J I— I (\j f-» (\4 (\i (\j 

OOOOOOOOOOOOOOOO 
+ ++ + + + 4-+-f-f+ + + + + + 

LU UJ UJ UJ UJ UU UJ UJ UJ UJ UJ LU UJ UJ UU UU 

in r-i O' ^ O' r-- <Nj ^ cNj CO CO cNj csJ ^ 
CO c\j (\j <M CO O" h* (\j in o >0 (\j >t 

CO r- cMn o r-- O' <N O' o o >o CM CM 

O' >r — « cn ^4“ >o nj o 0 ^ >r >o >3’ o m 

oj cn •-* (M ^ (\J CM *~i rsi (M 



f— icM fvj CM fSi fvJ nj rvj (M 

cncn >r ^4“ in in %o nO CO 00 C7" O' o o 

M M M Csl Kl »sj Csl M NJ M fSJ iSl NJ NJ K4 M 

3C3t J: Z JC 3C 'X3 3 2 3: 3E 2 ^ 
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COMPUTED DATA FOR SINGULAR LINES OF FIGURE 



a: 



OJ 


o 


o 


oooooooooooo 


oo 


o 


H- 


o 


o 


oooooooooooo 


oo 


o 


UJ 










•f 


s: 


LU 


UJ 


UJ UJ UJ UJ UJ UJ UJ UJ UJ UJ UJ UJ 


UJUJ 


UJ 


< 


o 


o 


oooooooooooo 


oo 


o 




o 


o 


oooooooooooo 


oo 


o 


< 


o 


o 


oooooooooooo 


oo 


r-i 


Q. 


o 


o 


oooooooooooo 


oo 


o 




o 


o 


oooooooooooo 


oo 


o 


C 

cc 


• 


• 
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